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We study the effects of doping a Mott insulator on the honeycomb lattice where spins interact 
via direction dependent Kitaev couplings Jk, and weak antiferromagnetic Heisenberg couplings J. 
This model is known to have a spin liquid ground state and may potentially be realized in correlated 
insulators with strong spin orbit coupling. The effect of hole doping is studied within a t-J-J^ model, 
treated using the SU(2) slave boson formalism, which correctly captures the parent spin liquid. We 
find superconductor ground states with spin triplet pairing that spontaneously break time reversal 
symmetry. Interestingly, the pairing is qualitatively different at low and high dopings, and undergoes 
a first order transition with doping. At high dopings, it is smoothly connected to a paired state 
of electrons propagating with the underlying free particle dispersion. However, at low dopings the 
dispersion is strongly influenced by the magnetic exchange, and is entirely different from the free 
particle band structure. Here the superconductivity is fully gapped and topological, analogous to 
spin polarized electrons with p x + ip y pairing. These results may be relevant to honeycomb lattice 
iridates such as A2lrC>3 (A = Li or Na) on doping. 

PACS numbers: 



I. INTRODUCTION 

The interplay of electron correlations and strong spin- 
orbit coupling (SOC) is currently attracting much atten- 
tion. Mott insulators with strong SOC, such as tran- 
sition metal oxides (TMO) of 5d elements, can display 
entirely different properties from those with weak SOC, 
such as cuprates, manganites and nickelates^l For exam- 
ple, the breakdown of the spin rotation symmetry allows 
for magnetic Hamiltonians very different from tradition- 
ally studied SU(2) symmetric models. This can introduce 
a new source of frustration^! leading to quantum spin liq- 
uid ground states. The Kitaev honeycomb lattice model, 
with spin dependent interactions between spin half mo- 
ments, is a remarkable example that admits an exact spin 
liquid ground stated It has recently been argued to be 
a natural Hamiltonian for a class of strong SOC mag- 
nets, such as the layered iridates A2&O3 (A = Na, LipEl, 
where Iridium atoms form the sites of a honeycomb lat- 
tice. In the iridium oxides, when an octahedral cage of 
oxygen atoms surrounds an Iridium ion, a j — 1/2 dou- 
blet is proposed on the Ir site^, for which a single band 
Hubbard model with strong spin orbit couplings can be 
invoked. In the Mott insulator, Ref. 17181 proposed that 
the spin couplings include both the isotropic Heisenberg 
term and the strongly anisotropic Kitaev coupling: 
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where SfS® is Ising coupling of the spin component a(= 
1, 2, 3) according to the type of (ij) boncP (see Fig. 1(a) |. 

Numerical calculation^^ f Eq. |l]) indicate the Kitaev 
spin liquid phase appearing at J = persists in the range 
< J < Jk/8. Although both A 2 Ir0 3 (A = Na, Li) are 



found to be magnetically ordered ^ 5 * 10 !, their transition 
temperatures are relatively low. Recent experimental pa- 
pers reporting magnetic susceptibility^^ have suggested 
that these iridates, particularly Li 2 Ir03 may be proxi- 
mate to the Kitaev spin liquid phase^LUJ Fits by exact 
diagonalization of the model Eq. ([I]) have reached similar 
conclusions^, but indicate that farther neighbor inter- 
actions also play a role. On the other hand, Ref. ll3llH 
proposed a rather different magnetic Hamiltonian, aris- 
ing from large trigonal distortions, and Ref. 15 proposed a 
quantum spin-Hall insulator. Future experiments should 
pin down the magnetic Hamiltonian in these materials. A 
different realization of the Hamiltonian Eq. Q is in per- 
ovskite iridate heterostructures of SrlrOg 16 , which pro- 
duces a honeycomb lattice when grown along the (111) 
direction. 

Motivated by these potential experimental realiza- 
tions, here we will study the effects of doping the 
Heisenberg-Kitaev model, and investigate the conduct- 
ing state that arises. To describe the physics of doping, 
we introduce the t-J-Jx model, with the hole doping of 
S per site, 
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where the projection operator V removes doubly occu- 
pied sites, and the chemical potential [i is adjusted such 
that c^Cia) = 1 — 5. The hopping term is nearest 
neighbor and spin independent. The symmetry of the 
honeycomb lattice along with reflection in the plane for- 
bid a spin dependence in the nearest neighbor hopping, 
as evidenced by microscopic considerations^. Farther 
neighbor hoppings can be spin dependent, but are ex- 
pected to be smaller and omitted in this minimal model. 
However, the spin-orbit interactions are nevertheless re- 



2 



tained in the Jk term. Similar Hamiltonian is also stud- 
ied in Ref.HUH 

The t- J- Jk model allows us the unique theoretical op- 
portunity of doping a magnet which is exactly soluble 
in the insulating limit (at the Kitaev point), and in a 
spin liquid phase. The exact solution singles out the cor- 
rect low energy variables — spins represented by neu- 
tral fermions (spinons), naturally motivating a slave bo- 
son formalism. Unlike in other studies of doped Mott 



insulator d 17 ! 18 !, here such a formalism can be a priori jus- 
tified. 

Our key results are as follows: (i) Doping Kitaev spin 
liquid leads to a spin triplet superconductor which spon- 
taneously breaks the time reversal symmetry, (ii) A first 
order transition occurs within the superconducting phase 
on increasing doping, which separates the two regimes 
SCi and SC2. In contrast, in a similar treatment of the 
well known square lattice t-J model, d-wave supercon- 
ductivity appears across the entire doping range at low 
temperature, and only quantitative properties are mod- 
ified with doping, (iii) In the low doping regime (SCi 
phase), quasiparticle dispersions are controlled by the 
magnetic exchange, and leads to a time-reversal-broken 
triplet superconductor with the same properties as a spin- 
polarized p x + ip y superconductor, which is fully gapped 
in the bulk but have chiral edge states and isolated Majo- 
rana modes in the vortex core^. This peculiar supercon- 
ducting state arises because of the unusual spinon disper- 
sion of the Kitaev spin liquid, (iv) At higher doping (SC2 
phase), the superconductor obtained reflects the bare dis- 
persion of electrons, and can be smoothly connected to 
the weak coupling limit, where magnetic interactions lead 
to pairing near the Fermi surface. 

This paper is structured as follows. We begin by an- 
alyzing the quantum order underlying the Kitaev spin 
liquid, characterized by the symmetry transformations 
of fractionalized excitations, a description known as the 
projective symmetry group^l (PSG). We find that the 
Kitaev quantum order locks the spin and gauge rota- 
tions together; the two holon species transform like a 
spin, and spontaneously break time reversal when con- 
densed. Next we map out the mean field phase diagram 
within the SU(2) slave boson formalism as constrained 
by the Kitaev PSG, exact at zero doping, and demon- 
strate that the SCi and SC2 phases are dominated by 
different physics. Controlled by the quantum order, a 
time-reversal-broken triplet superconductor SCi emerges 
from the doped Kitaev spin l iquid . We close with com- 
ments on related recent worlipLLM] anc j ^he relevance of 
our result to experimental realizations. 



would like to analyze the symmetry property of the model 
and its spin liquid ground state. 



A. Symmetries of the Kitaev Model 

First, we consider the space group symmetries of the 
model. The symmetries are most naturally expressed by 
embedding the honeycomb within a 3D cubic lattice, ex- 
actly in the same manner that the Kitaev honeycomb 
model arises in three-dimensional layered iridates. Then 
the symmetry transformations, which due to spin-orbit 
coupling act simultaneously on spin and space, are rep- 
resented in the same manner on the spin space and the 
3D real space. 

Specifically, the space group is generated by two trans- 
lations Ti and T 2 , an operation C 6 composed of a 6-fold 
c-axis rotation followed by a reflection across the lattice 
plane (the c = plane), and a reflection a across the 
x = y plane, as illustrated in Fig.[T] 
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(b) 



FIG. 1: (Color online.) (a) Lattice symmetries of the Kitaev 
model. The operations Ce and a act simultaneously on lat- 
tice and spin. The three bond types (a = 1,2,3) are colored 
red, green and blue respectively, (b) A hexagon plaquette 
embedded in the cubic lattice. The c-axis is the (111) axis. 
The 6-fold c-axis rotation is not a symmetry by itself, but be- 
comes a symmetry when combined with the reflection across 
the lattice plane. 



Besides the space group symmetries illustrated above, 
the Kitaev model is also symmetric under time reversal 
T ■ Time reversal has no effect on the lattice but acts 
as icr 2 followed by complex conjugation K. on the spins. 
While T 2 = — 1 on a single spin, the global time reversal 
symmetry operation acting on the bipartite honeycomb 
lattice squares to +1. Combining T with the space group 
yields the full symmetry group (SG), with the presenta- 
tion SG - (r.Ti.Ta.CB.trlT 2 = I, a 2 = 1, (C 6 ) 6 = l) 
subject to 13 definition relations, listed in Eq. (A9). 



II. KITAEV SPIN LIQUID 

To explore the physics of the t-J-Jji model, we start 
from the well-controlled undoped and J = limit, where 
the model reduces to the Kitaev model. Its exact solution 
is given by Kitaev^ and is already well-known. Here we 



B. Symmetries in a Schwinger Fermion 
Decomposition: the Projective Symmetry Group 

In order to study the Kitaev spin liquid and nearby 
phases, we must decompose the spin operator S" — 
\fi°afi into fermionic spinons fj = (Z^,/^), with a a 



3 



being the Pauli matrices. Compared to the spin oper- 
ators S", the spinon operators f%& have an additional 
SU(2) gauge structure, best seen by arranging the oper- 
ators into the following matrix-^ 



Ft 




(3) 



Any right SU(2) rotation F % -> FiG : G € SU(2) leaves 
the physical spin Sf (and hence the spin Hamiltonian) 
unchanged, as can be seen from the following equivalent 
expression of Sf 



1. 



Sf - jTrf V,,. 



(4) 



Therefore the right rotation G corresponds to a gauge 
SU(2) rotation, whose generators (the SU(2) gauge 
charges of spinons) are given by 



(5) 



On the other hand, the left rotation F^ WFj : U € 
SU(2) corresponds to the spin SU(2) rotation, whose gen- 
erators are the spin operators S". 

Because of the gauge SU(2) redundancy in the 
Schwinger fermion representation, any SU(2) gauge op- 
eration leaves the physical spin system invariant. Any 
operator acting on the spins, such as a symmetry trans- 
formation, may also act within this SU(2) gauge space. 
Thus when we fractionalize spins in a Schwinger fermion 
decomposition, we must also specify how the symme- 
try operations of the model act within the gauge free- 
dom. This extra information, known as the projective 
symmetry grourJ^ (PSG), characterizes the fractional- 
ized phase. Symmetry operations therefore consist of a 
symmetry group operation g £ESG with the correspond- 
ing spin operation U g and gauge operation G g , such that 
the spinons transform as 



F^Ul^F^G^i). 



(G) 



The index i labels the site. 

In fact, the spin operation U g (i) 



U„ are always site- 



independent, so the site index may be omitted. f/ fl 's are 



given by 



U Tl = U T2 = 1, 

UcM) = Uc e (B) = <?c 6 
U a (A) = Ua(B)=a a , 



(7) 



where ac e = (co + i&i + i&2 + icr$)/2 and a a = i{<j\ — 
o'2)/v / 2- Co is the 2x2 identity matrix. These matrix 
representations are literally translated from the descrip- 
tions of the symmetry operations on the cubic lattice, see 
Fig. 1(b) The antiunitary time reversal operation can 



be represented by a unitary transformation followed by 



a complex conjugation /C, which transforms the spinons 
by 

Ft -> ICU^FiGr^lC, (8) 

where the unitary operation acting on the spin reads 

U T {A) =U T {B) = ia 2 . (9) 

The complex conjugate operation K, flips the sign of the 
imaginary unit, i.e. JCi = — ifC, while keeping everything 
else invariant (K? = 1). 



C. Projective Construction for the Kitaev Spin 
Liquid 

The Kitaev model can be solved exactly^ by intro- 
ducing 4 Majorana fermions xf ( a — 0, 1, 2, 3) on each 
site, and rewriting the spin operators as Sf = iXiXf 
under the constraint XiXjXiXi = 1/4- The Majorana 
fermions are normalized as |y?,x" } — ^w^aa 1 in this 
work. It has been pointed oupSthat under certain SU(2) 
gauge choice, the Majorana fermions xf are related to the 
Schwinger fermions fi a by the following matrix identity 
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or more explicitly as 

,21 



XD- The Majorana fermions introduced by Kitaev are 
just another representation of the spinons. All the emer- 
gent SU(2) gauge structure for Schwinger fermions fia- 
applies to the Majorana fermions xf as well. 

The exact ground state can be obtained by the follow- 
ing projective construction^. First take the Majorana 
bilinear Hamiltonian 



H = J K J2«jXiX°j 

(ij) 
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where the bond parameters ufj = (ixfxf) i a — 0, 1, 2, 3) 
can be regarded as the mean field ansatz, self-consistently 
given by 



-0.262433 if a = 0, 
1/2 ifa = o, 

otherwise. 



(12) 



Here a denotes the type of the bond (ij). We choose 
i E A sublattice and j £ B sublattice to be the posi- 
tive bond direction. Given the ansatz Eq. ( |12[ ) , the mean 
field Hamiltonian Eq. ( |Tlj ) produces a graphene-like band 
structure for x° and degen erate flat bands for x ■> X 2 an< l 



X , as shown in Fig. 2(a) Take the Majorana Fermi liq- 



uid ground state and project to the physical Hilbert space 
by imposing the condition xlXiXiXi =1/4, the resulting 
state is the exact ground state given by Kitaev. 

The spin correlation in this state was shown to be 
short-rangecp3, which identifies the ground state of the 
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(a) 



(b) 



FIG. 2: Mean field band structure of Majorana spinons (a) 
in the undoped limit, (b) with doping St/ Jk = 0.2. The inset 
shows the x500 zoom-in around the K point. 



Kitaev model as a quantum spin liquid. However what re- 
ally differentiates the spin liquid from a trivial spin disor- 
dered paramagnetic state is the quantum order^ encoded 
in the Majorana Fermi liquid from which the spin liquid is 
obtained by projection. Given the particular mean-field 
ansatz parameterized by ufj, the \° fermion has a band 
structure different from x 1,2 ' 3 j so it is no longer possible 
to mix x° with the other Majorana fermions. Thus the 
emergent SU(2) gauge structure of mixing spinon flavors 
is broken down to the Z2 gauge structure of changing sign 
of x a - The broken gauge structure can be imagined as 
a hidden order of spinon superconductivity^!. Although 
it will not manifest as electron superconductivity in the 
spin liquid due to the lack of charge fluctuation, its ex- 
istence as a quantum order is real, and will be revealed, 
once the charge fluctuation is introduced by doping. 



D. Projective Symmetry Group of the Kitaev Spin 
Liquid 

More precisely, the quantum ordei^ of the Z 2 spin liq- 
uid is characterized by the PSG of the mean field ansatz. 
The PSG of Kitaev spin liquid can be determined start- 
ing from the fact that x° is a special flavor which should 
not be mixed with other flavors, any PSG operation must 
at least preserve the flavor of x°- X° appears in the F 
matrix as F ~ X 0(J o, while F transforms under PSG op- 
erations as F — > W g FGg, so apart from some sign factor, 
X°co ^X^U^Gg. Therefore, to preserve the flavor of 
X°, G g = ±U g is simply required to hold for all g £ SG: 
the gauge operation G g must always follow the spin op- 
eration U g up to a sign factor. From the spin operations 
U g given in Eq. Q and Eq. |9| it is not difficult to figure 
out the gauge operations G g , which read 



Gti — Gt 2 — 1, 
G C6 (A) = -G Ce (B) = <?c 6 , 
G a {A) = ~G a {B) = a ai 
G T {A) = -G T (B)=i<r 2 . 



(13) 



The matrices ac 6 and a a were defined right below Eq. ^ . 
The sublattice-dependent sign factors are determined as 
follows. Both Cq and a switch the sublattice A and 



B, carrying u\ B to u% A under the lattice transforma- 
tion. However, u% = —U% is odd under the reversal 
of bond direction, so in order to keep it unchanged, 
the sign must be rectified by the gauge operation that 
follows, therefore both Gc 6 and G a have a sign differ- 
ence between the sublattices. However for the time re- 
versal operation, under complex conjugate i — > —i, so 
u ab = (*XaXb> HxaXb) = ~u?j, thus the gauge 
transform Gj- must also carry the sublattice-dependent 
sign to compensate the sign generated by the complex 
conjugate. 

A prominent property of the PSG of the Kitaev spin 
liquid is that U g and G g are always the same (up to 
a sign), which implies that the spin and gauge degrees 
of freedom are locked together by the underlying quan- 
tum order in the spin liquid state. As a result, the PSG 
operation U^FG g literally carries out the rotations and 
reflections by treating \° as a scalar and x = (x 1 ) X 2 ■> X 3 ) 
as a pseudo vector. Therefore Cg actually permutes 
X 3 - ^ X 2 X 1 ~~ * X 3 j an d o exchanges x 1 X 2 > with 
some additional sign factors (see Tab.|l]), thus giving ex- 
actly the right transforms to preserve all the mean field 
ansatz, which can be checked straightforwardly. 

TABLE I: The PSG transforms of Majorana fermions. 



Xa 
XX 

x\ 

X 3 A 
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T 



XT 
X% 

xl 



xT 

Xa 

XA 
XA 



1^ 

xl 

Xb 

xl 



xT^xT 

Xs -Xa 
X% -Xa 
Xb -Xa. 



~^b~ 

-xl 
-xl 
-xl 



-xJ 

XA 
XA 

xl 



Xa 
Xa 
Xa 



-xJ 
-xl 
-xl 
-xl 



In conclusion, the PSG of the Kitaev spin liquid is de- 
fined by Eq. (|6| in general (and by Eq. Q for the time 
reversal operation), with the spin and gauge transforms 
specified by Eq. 0, Eq. ^ and Eq. Q. Its effect on 
the Majorana spinons are concluded in Tab.|lj This PSG 
belongs to the class (I)(B) according to the PSG classi- 
fication of Z2 spin liquid on the honeycomb lattice (see 
Appedix |A| for details of the classification). 

All the PSG's in this class have the common property 
that the gauge charge is reversed under time reversal just 



the same as the spin. To see this, substitute Eq. ( 10 ) into 
Eq. Q and Eq. and write the spin and gauge charge 
operators in terms of Majorana fermions as 



; (Xi Xi 



Ki = -(x iXi 



Xi x Xt 



(14) 



:Xi x Xt 



Applying the PSG transformation rules of the time re- 
versal: xl ~^ X% Xb ~Xb ( see Tab.|l| and i — > —i, 
it is easy to show that both the spin and gauge charge 
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operators are odd under time reversal 

Si^-SuKi^-Ki. (15) 

Therefore, there are in principle two ways to to break 
the time reversal symmetry in the Kitaev spin liquid: 
one is to polarize the spin and the other is to condense 
the gauge charge. The spin polarization can be achieved 
by applying an external magnetic field in the (111) direc- 
tion (perpendicular to the lattice plane), which drives the 
gapless Kitaev spin liquid into the gapped non-Abelian 
phase^l. In the following, we will explore the second 
possibility, namely the gauge charge condensation. This 
can be achieved by introducing the gauge charge through 
doping the spin liquid. According to the SU(2) slave bo- 
son theory, the condensed holon will pick out an SU(2) 
gauge direction and break the time reversal symmetry 
spontaneously. 



Let 1 0) s i a vo be the vacuum state of both spinons and 
holons, s.t. / i(T |0) s i avo = 6i„|0) s i avc = 0. Then on each 
site, there are only three physical states in the Hilbert 
space: 

|0) = ^(4+4/i/i) I0)slave, 

c t r |o> = /t r |o) slave , ( 18 ) 

i|0>=/i|0) s lavc. 

Here |0) denotes the electron empty state. The double 
occupied state is automatically ruled out from the phys- 
ical Hilbert space in the SU(2) slave boson formalism. 

Each empty site has one holon, therefore the doping 8 
is: 5 = -h + b\ 2 b i2 ), where N denotes the total 

number of sites. Adopting Gutzwiller approximation, the 
spin operator will be written as Sf — iXiXifi — £)• 



III. DOPING THE KITAEV MODEL WITHIN 
SU(2) SLAVE BOSON THEORY 

A. SU(2) Slave Boson / Schwinger Fermion 
Representation 

We now consider doping (say) holes into the insulating 
magnet, while preserving the strong onsite correlations 
that penalize double occupancy. As discussed above, the 
exact solution of the Kitaev spin liquid is naturally ex- 
pressed within a particular kind of Schwinger fermion / 
slave boson representation. The most naive way is to 
directly assign the spinons to the electrons Ci a = bifi a 
with additional U(l) slave boson bi to carry the elec- 
tric charge. However this approach completely neglects 
the SU(2) gauge redundancy in the spin liquid: annihila- 
tion of a spin up electron by can be accomplished (in 
the spin sector) either by the annihilation of up spinon 
f-f or by the creation of down spinon /j (to neutralize 
the up spin into spin singlet), so t he ele ctron operator 
must be a linear combination of bot h 1 18 * 30 -!, formulated as 

Cit = ^$i/*t - *l = jsOtiM + 4/1) 

equivalently as^ 



or 



V2 



(16) 



where Cj, Fi and Bi are 2x2 matrices of operators 



j 



, Bi 



4 
4 



-ba 
hi 



(17) 



and Fi is given by Eq. ^ in terms of Schwinger 
fermions or equivalently by Eq. ( 10 ) in terms of Majo- 
rana fermions. The holon creation operators b\ x and b\ 2 
carry different SU(2) gauge charges, but the same electric 
charge as a hole Cj. 



B. SU(2) Gauge Charge 

As both spinons and holons carry the SU(2) gauge 
charges, the gauge SU(2) generators K\ (I — 1,2,3) are 
generalized from Eq. ((1) to 



K 



1. 



1. 



Tr FicriF- - -Tra 3 B <nB 



(19) 



or explicitly written as (with implicit sum over dummy 
indices) 

K\ = -\{iX°iX\ + \e lmn xTx n i + b lv a l vv ,b\ v ,\ (20) 

where ei mn is the Levi-Civita symbol. It can be veri- 
fied that \K\ , F^ = \F i( j u [K\,Bi] = -\oiBi, therefore 
[Kid] = 0, showing that K\ are indeed the generators 
of gauge SU(2) transforms that leave the electron opera- 
tors unchanged. 



The physical state, as enumerated in Eq. (18), are 
SU(2) gauge invariant. Therefore the SU(2) singlet con- 
dition K\ — should be imposed. This condition is 
equivalent to the single occupancy condition for both 
spinons and holons, as is evidenced from Kf = (1 — 

/1/it - &fu ~ b\ x b lX + b\ 2 b l2 )/2 = 0. 

The PSG operations are naturally extended to the 
holons, such that they transform as 



Bi 
Bi 



B 



ff(»)' 



(for g ? T) 



KGiyBiK. 



In particular, under the time reversal operation, 



(21) 



(22) 



One can see the holon SU(2) gauge charges transform 
under time reversal in a way similar to the physical 
spins. Therefore one could expect that the condensa- 
tion of holons will spontaneously breaks the time reversal 
symmetry. 



( b A1 ^ 




~bA2 \ 


I b A2 i 
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C. Mean Field Phase Diagram 

The exact solution of the Kitaev spin liquid at zero 
doping involves an enlarged hilbert space with spinons 
and holons which implements a particular PSG. We ex- 
pect these deconhned excitations, which transform under 
symmetry operations as defined by the Kitaev-limit PSG, 
to survive into finite doping. At small finite doping the 
SU(2) slave boson mean field with this particular PSG 
becomes inexact, but should still provide the most accu- 
rate treatment possible. 

Using Eq. (16), the £-J-Jr model can be written in 
terms of spinons and holons (for simplicity we set J = 0, 
finite J is discussed in Appendix [C]) . Then use the mean 
field treatment by introducing the following mean field 
parameters: 



(23) 



we arrive at the mean field Hamiltonian (see Appendix 
[Blfor detailed deductions) 



Hmf = £ U^xfx" + Wfriblbj* + h.c.) 



(ij) 



(24) 



where summation is implied over repeated indices a = 
Q . . . 3, v = 1, 2 and Z = 1 ... 3. The hopping amplitudes 
for fcrmions and for bosons WK should be determined 



self-consistcntly from 



U a 



t 2 

4ZZ« + c - c -) 



+ J K (l - S) 2 (u^S 0a + tljioa), 



(25) 



a=0 



The index a denotes the direction of (ij). The boson 
chemical potential fi is chosen such that J2i v^iv ™) = 
6N '. The SU(2) gauge charge operators are given in 
Eq. (20 1. The gauge potentials a \ are chosen to enforce 
the SU(2) gauge singlet constraint on average (K\) = 0. 
In the undoped limit, Eq. (24) reduces to the mean field 
description of the spin liquid exact solution. With fi- 
nite doping, the hidden superconductivity of spinons will 
be rendered into the true superconductivity of electrons 
once the holons condense. 

We would like to stress that the quantum order of the 
Kitaev spin liquid puts a strong constraint on the possible 
form of the mean field ansatz. This quantum order is 
described by the Kitaev spin liquid PSG as discussed 
previously. We assume that this PSG is respected by 
the mean field solution throughout, and that symmetry 
breaking occurs only through holon condensation. At 
small dopings this is required by continuity to the Kitaev 



solution. The most general parameterization of the mean 
field ansatz under the PSG restriction is as follows. First 
assign on the type-3 bond 



w 



Uq,U. 



Ub, U; 



(26) 



Then the mean field parameters on the other bonds are 
obtained by using the PSG operation to carry the above 
assignment throughout the lattice, uq, u a , and w are 
real numbers that parameterize the mean field ansatz. 

Based on the parameterization, a self-consistent mean 
field solution of Eq. ( 24 ) gives the phase diagram shown 
in Fig. [3] We show results for J = 0. Introducing J < 
Jk/8. to remain within the boundary of the spin liquid 
phased, has little effect on the phase diagram (discussed 
in Appendix [Cl. 




St/J K 
Chern 
number 



FIG. 3: (Color online) Mean field phase diagram for t = IOJk 
and J = 0. The low doping Kitaev spin liquid (KL) phase 
and the high doping Fermi liquid (FL) phase are separated 
by first order transition. Once holons condense, two classes 
of superconducting (SCi and SC2) phases appear. The bar 
below shows the Chern number of the superconducting state. 



D. Spin Liquid and Adjacent Phases 

In the undoped limit, one recovers the Kitaev spin liq- 
uid mean field parameters: Ub = w = 0, and uq and u a 
arc determined by the following self-consistent equations 



1 /3JkUq 
- tann . 

2 2 



u a = 



± V |F( fc )|tanh^M. 



(27) 



where T(k) = e ifc » + 2 e - lfe «/ 2 cos^W 2 ), and N is the 
number of sites. At zero temperature, the solution is 
u = —0.262433 and u a = 1/2, corresponding to the 
exact ground state of the Kitaev model. So the SU(2) 
slave boson mean field theory is asymptotically exact in 
the small doping limit. At the mean field level, a finite 
temperature transition is found at T c — Jk/^, above 
which [T > T c ) all the mean field parameters vanish, 
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Uo = u a = Ub = w = 0. The confining gauge fluctu- 
ation will recombine spinons and holons into electrons, 
resulting in a paramagnetic (PM) phase. 

With increasing doping, mean field parameters uj, and 
w grow in proportional to 8, and eventually trigger a first 
order phase transition at 8 C ~ 2u a J\^/t, see Fig.g The 
transition is driven by the competition between the ki- 
netic energy of holes (t term) and the magnetic energy of 
spins (Jk term). The magnetic energy favors the Kitaev 
spin liquid state, in which the mobility of x 1,2 ' 3 fermions 
is sacrificed (as they form degenerate flat bands). For 
larger doping, more kinetic energy can be gained by al- 
lowing ^ 1 ' 2 - 3 fermions to move in the same way as x , 
as uo ~ u a ~ U b, so that the flat band gets dispersed 
as shown in Fig. |2(b)| In the large doping limit, all fla- 
vors of Majorana fermions move with the same ampli- 
tude, providing identical graphene-like band structures, 
which can be recombined into band electrons, labeled by 
Fermi liquid (FL) in Fig. [3] As discussed below, the na- 
ture of superconductivity is very different depending on 
the normal state, Kitaev spin liquid or FL, from which 
it emerges. 



u 

0.6 fKL 
0.4 
0.2 
0.0 
-0.2 



FL 



0.( 1.0 1.5 2.0 2.5 3.0 



St/J K 



where n = (-y/3/2, -1/2), r 2 = (^3/2, -1/2), r 3 = 
(0,-1) denote the three displacement vectors from site 
A to site B, and b labels the singlet (b — 0) or triplet 
(b = 1,2,3) channels. o>b — z J e<JbZ with z = (zi,Z2) J 
refer to the holon condensate amplitude, and d a b — 
uq — u a + 2(u a — Ub)5 a b parameterize the the spinon pair- 
ing amplitude. The electron superconductivity is a joint 
effect of holon condensation and spinon pairing. 

Obviously aa — for whatever z, so A A b,o — 0, thus 
the electron paring is purely triplet. This demonstrates 
the spin-gauge locking effect of the Kitaev spin liquid, 
that a singlet in the spin space will be rendered by the 
PSG to a singlet in the SU(2) gauge space (seen from the 
expression of a&). However gauge charges can not con- 
densed to a singlet state due to their bosonic nature, thus 
the single pairing is ruled out, as long as the quantum 
order persists. 

The superconductivity transition temperature in the 
phase diagram shown with the dashed line is estimated 
as follows. At small doping, the phase stiffness pb = tbS 
is proportional to doping, where tb = M(uq + u a + 
2ub)/8, and T c is estimated from the Kosterlitz-Thouless 
transition^ temperature T c = Trpb/2. At large doping, 
the mean field gap is small which controls T c ~ Af, where 
Af ~ J K (uo + u a) 1/2 (l - <5 2 )/ 4 - In between, we inter- 
polate via the formulcPT" 1 = (irp b /2)- x + Aj 1 . Note, 
due to the absence of a finite temperature transition of 
the two dimensional free bosons, a naive mean field tran- 
sition temperature is not specified. 



FIG. 4: (Color online.) Mean field parameters uo, u a , Ub v.s. 
doping 8 at zero temperature. The arrow indicates the 1st 
order transition between KL/SCi phase and FL/SC2 phase. 
The calculation is done at t — IOJk and J = 0. 



E. Holon Condensation and Superconductivity 

At low temperature, the holons condense to their band 
minimum at zero-momentum, leading to the following 
condensate amplitude (y = 1,2) 



{b Av ) = z v , (b BlJ ) = iz v , 



(28) 



with the density \z\\ 2 + \z2\ 2 = S following the doping 
level. The electron pairing is found between opposite 
sublattices (because the intra-sublattice coupling of x" 
is forbidden by PSG): A A B,b{k) = c J kA eabC- k B, {b = 
0,1,2,3) where c kA{B ) = (cfe^(s)f, c fc A(B)i) T denote the 
electron operators in the momentum spa ce a nd e = itj2 is 
the anti-symmetric matrbP^l Using Eq. (16 1, the pairing 



is expressed in terms of the mean field parameters: 



(fe) 



3 

y 2^ dabe 

a=l 



(29) 



F. Symmetry and Topological Properties 

The mean field Hamiltonian of the Kitaev spin liq- 
uid appears surprising at first, since the only Majorana 
fermion with extended hopping is % , the real part of 
/♦, which seems to single out one spin species and break 
the time reversal symmetry. Actually, this is a gauge ar- 
tifact. The SU(2) rotations between and Ci a will re- 
store the time reversal symmetry on the electron level for 
the spin liquid. However, the SU(2) gauge redundancy is 
parameterized by holon fields bi V and must be resolved 
as the holon condenses. So, as has been discussed from 
the PSG prospective, the holon condensation must break 
the time reversal symmetry spontaneously, leading to a 
class D superconductoi!^!, denoted as SCi, with uniform 
magnetization (S) ~ z'az. 

Let us elaborate on the microscopic mechanism which 
gaps the x° Majoranas in SCi. If we view the charge and 
spin as separate excitations, one may expect the same 
spectrum as the Kitaev spin liquid, i.e. gapless x° Majo- 
rana modes, to persist into the superconductor. However, 
the time reversal symmetry, which protects this gapless- 
ness in the spin liquid, is lost in the superconductor. This 



can lead to an energy gap for x° ( a s shown in Fig. 2(b) I, 
tied to the strength of the condensate. Because the uni- 
form SU(2) gauge charge provided by the holon conden- 
sate offsets the SU(2) gauge potential a 1 (I — 1, 2, 3) from 



zero, in order to preserve the overall gauge singlet condi- 
tion. It is found that a 1 ~ SJk increases with doping. For 
small doping 6, we treat a'/Jx as a perturbation. Inte- 
grating out the gapped Majorana modesv 1 ' 2,3 generates 
next nearest neighboring (nnn) (Fig. 5(b) ) hopping of x° 
fermions through a 3rd order perturbation correction, as 
illustrated in Fig. 5(a) The effective Hamiltonian for x° 
reads 



(ij) 



(30) 



where v = a 1 a 2 a 3 /(8 J^ u q) an d ((ij)) denotes the ori- 
ented nnn bond, with the bond direction specified in 



Fig. 5(b) According to the Kitaev spin liquid PSG (see 
Tab.|l|, the nnn coupling term is time reversal odd (since 
i — > —i), and is allowed only because time reversal sym- 
metry is broken by the gauge charge condensation here. 

The resulting x° Hamiltonian Eq. ( f30| ) is a Majorana 
version of the Haldane model 34 . It is known that the 
nnn coupling gaps the Dirac cones and leaves one unit of 
Chern number in the ground state. This requires all a 1 
to be nonvanishing. It is actually energetically favorable 
for the holon condensate (i.e. magnetization) to be in 
the (111) (or equivalent (±1 ± 1 ± 1)) direction, corre- 
sponding to a 1 = a 2 = a 3 which maximizes the spinon 
gap m = Therefore in the small dop- 

ing limit, the ground state is a fully-gapped topological 
superconductor with +1 Chern number, which implies a 
gapless chiral Majorana edge mode and a Majorana zero 
mode in the vortex core. This is the same topology as 
a p x + ipy superconductor of spin polarized fermions^ 
here the "spin-polarization" arises from the peculiar dis- 
persion of fermions in the Kitaev spin liquid. At larger 
doping the Chern number changes, as shown in Fig. [3] 
The transition +1 — » —2 in the SCi phase corresponds 
to a band gap closing at M point due to the softening of 
X 1,2 ' 3 modes. 
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FIG. 5: (Color online.) (a) The path of 3rd order pertur- 
bation. The 4 Majorana fermions on each site are denoted 
by their flavor indices. The effective second nearest neighbor 
hopping of x° fermion is bridged by two nearest-bond hopping 
of X 1 ' 2 ' 3 - The on-site flavor changing process is assisted by 
the time-reversal-broken gauge potential, (b) Gray dashed ar- 
rows indicate the directions of the second nearest neighboring 
bond. 



IV. DISCUSSION AND CONCLUSION 

A. Overdoped Regime and Weak Coupling BCS 

In the overdoped FL phase where correlations are 
weak, the superconductivity (SC2) can be studied un- 
der the BCS paradigm by treating H aK as an interac- 
tion and decomposing it into the Cooper channel. In 
the small J limit, The instability is found in the spin- 
triplet pairing channel, because the spin model is ferro- 
magnetic. To first order in weak coupling both the time 
reversal invariant superconductor (the two dimensional 
analog of He3 B phase) and the time reversal symmetry 
broken triplet superconductor (the analog of the He3 A 
phase) are degenerate. To next order, the calculation 
in Ref J21 showed that the time-reversal-invariant p-wave 
superconductor is preferred. Beyond weak coupling it is 
hard to decide which of these two possibilities is realized, 
a problem that is well known from He3 physics^. Here, 
our choice of PSG selects the time reversal (T) broken 
state, while a different choice would yield the T symmet- 
ric state. Therefore we mention both these possibilities 
as potentially relevant to the material at hand at high 
doping. In either case, the SC2 phase is dominated by 
the Fermi liquid physics and is separated by a first order 
transition from the spin-liquid-controlled time-reversal- 
broken SCi phase elaborated in this work. Because of 
the distinct underlying mechanism, its is not surprising 
that SCi and SC2 can be quite different in many aspects. 



B. Conclusion 

A time-reversal-broken spin-triplet topological super- 
conductor was found in the doped Kitaev spin liquid 
within the SU(2) slave boson formalism. A first or- 
der quantum transition around S c ~ JkA separates the 
spin triplet superconductor into two distinct classes: SCi 
(controlled by Jk) is governed by the spin liquid physics 
and reflects the underlying quantum order, while SC2 
(controlled by t) is a more conventional BCS-type su- 
perconductor. Although both ultimately trace their ori- 
gins to the magnetic couplings, the detailed mechanisms 
are rather different. This is in sharp contrast to the t-J 
model in the context of cuprates, where, at least qualita- 
tively, d-wecve superconductivity is realized throughout. 

A pro mising candidate material is A2lr03 (A = Na, 
Li) 4 * 5 * 10 *, although experiments suggest magnetic ground 
state, rather than spin liquid. However, it has been ar- 
gued that doped charges are more mobile in spin liquids, 
as compared to antiferromagnetic states where they in- 
terfere with the ordered patter rP. Therefore one may 
hope that the results derived here also hold for magnetic 
ground states that are proximate to the Kitaev phase. 
Our main prediction is that doping these systems should 
lead to spin triplet topological superconductors with su- 
perconducting T c a fract ion of the magnetic exchange. 
Assuming Jk ~ 100- 150KP^E^a crude estimate of maxi- 
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mum superconducting transition temperature is 15- 20K. 
Although we are not aware of doping studies on this class 
of materials, t he re lated iridium perovskite Sr 2 Ir04, a 5d 
cuprate analog^SI has been doped in the buUP^l, and re- 
cent years have witnessed significant progress in doping 
techniques. We hope our results will spur future experi- 
ments in this direction. 



T2, Cq and a with the following 13 definition relations 
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1, 



TTiTTf 1 = TT 2 TT^ = 1, 

CrTiCT- 1 ^ 1 = oTi<J- X T^ 1 = 1, 

T 2 = C% = a 2 = 1, 
TCeTCg- 1 = 1, 
7V7V- 1 = 1, 
C§oC§o = 1. 



(A2) 
(A3) 
(A4) 
(A5) 
(A6) 
(A7) 
(A8) 
(A9) 



In general each definition relation takes the form of 
• • • 5251 = 1, where • • • g 2 gi denotes a sequence of symme- 
try group operations. Then according to Eq. |6]), under 
the PSG operation, the spinon matrix Fi transforms as 



F > "> u l u l ■ ■ ■ F - g *9i(i) ■ ■ ■ G 92 (ffi(0) G gi (i). 



(A10) 



Because the bunch of operations • • • 3231 actually result 
in the identity operation, so they must not affect the spin 
degree of freedom: l/J U^ 2 ■ ■ ■ — <jq and must also restore 
the original lattice site: • • ■ g 2 gi(i) — h hence the PSG 
operation becomes a pure gauge operation 



p. '"JIM, 1 p. 



■■G B2 ( 9l (z))G gi (i). 



(All) 



Appendix A: Z2 Projective Symmetry Group on 
Honeycomb Lattice 



Here we present the classification of Z2 projective sym- 
metry group (PSG) on Honeycomb lattice without spin 
rotational symmetry (but preserving time several sym- 
metry). 144 solutions of algebraic PSG were found. 

On the Honeycomb lattice, each unit cell is labeled by 
its integer coordinates X\ and x 2 along the translation 
axes of T\ and T 2 . A spin site is further specified by its 
sublattice label A or B within the unit cell, see Fig 
The symmetry group operators act on the lattice by 



la 



Ti{x ll x 2 ) = {xi + l,x 2 ), 
T 2 (x 1 ,x 2 ) = {xi,x 2 + 1), 
C 6 (xi,x 2 , A) = (xi - x 2 .xi,B), 
C 6 (xi,x 2 ,B) = (xi - x 2 - l,Xt,A), 
a(x ll x 2l A) = (x 2 ,xi,B), 
a(xi,x 2 ,B) = (x 2 ,xi,A). 



The sublattice label is omitted if a formula holds in both 
sublattices. Later we will refer to the principal unit cell 
by omitting the unit cell index, i.e. (0, 0, A) = (A), 
(0,0, B) = (B). The representation of symmetry opera- 
tors in the spin space will be given after further discussion 
in Eq. 

The symmetry group of a general spin model on the 
Honeycomb lattice is generated by 5 generators T, T\, 



All the pure gauge operations that leaves the mean field 
ansatz invariant constitute a subgroup of the PSG, known 
as the invariant gauge group (IGG). So we must have 

• • • G 32 (gi(i)) G gi (i) e IGG. Here we are interested in 
the classification of Z2 spin liquid, so we will focus on 
the case that IGG = Z 2 . Thus for each definition relation 

• • • 3251 = 1, there is a corresponding PSG representation 



G g2 (g 1 {i))G g 



(i) = Vn 



(A12) 



where r) m = ±00 will be used to denote the sign factors 
hereon. For the 13 definition relations, we introduce 13 
sign factors r)x , 772, ■ • • , r]i3 to denote the corresponding 
IGG elements. In the following, we may write the PSG 
representation Eq. ( A12) in short as • • • G 92 G gi = rj m by 
omitting the site labels so as to save the space. 

However special attention should be paid to the time 
reversal operation, because it involves the complex con- 
(Al) jugate operator K, which does not commute with G g in 
general. As can be seen from Eq. (|8|), K. must be placed 
right after each G- 
be represented as 



right after each Gr- For example, TCqTC 6 1 = 1 should 



G r (*)/CG C6 (z')Gr(*')^G C6 1 (z') = mi , 



(A13) 



where i' — C 6 _1 (i). Here we have used the rule that 
G g -i(g{i)) = to simplify the inverse operations. 

To classify the PSG's one should take care of the gauge 
redundancy in the solution of G g . Two PSG's are gauge 
equivalent if their solutions of G g are related by a set 
of local SU(2) gauge transform G g {i) -> H / J ( . ) G g (i)W i : 
Wi G SU(2). To reduce the gauge redundancy, gauge 
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fixing will be used while solving the equations of G g . 
First of all, the relative gauge between the unit cells can 
be fixed by setti ng Gt 2 (xi,X2) = &o, and Gji(xi,0) = 
(Jo, then Eq. ( A2 1 can be represented as Gtj (xi, £2 + 1) = 
t]iGt 1 {xi,X2), which gives the solution for translations 



G Tl (x 1 ,x 2 ) = Vi\G T2 {x 1 ,x 2 ) = <t . 



(A14) 



Substitute Eq. ( jA14 1 into the PSG representation of 
Eq.dAl: G T lCG Tl lC = r) 2 G Tl and G T 1CG T2 G T 1C = 



V3GT2 j we obtain 



Gr{x\ + 1,X2)K,Gt{x\,X2)K = T]2, 

G T (x 1 ,x 2 + l)K,G T (x 1 ,x 2 )IC = 7/3, 



(A15) 



while on the other hand, from GtICGt IC — ? 7s, we know 
K.Gf{xi,x%)K = rjsG^- 1 (xi, X2), so Eq. (A15) becomes 



G T (xi + l,x 2 ) = rj2TigGr{xi,X2), 
G r (xi,x 2 + 1) = V3VsGr(x 1 ,X2). 



The solution is 

G T (x 1 ,x 2 ) 



(A16) 



(A17) 



Similarly by inserting Eq. (A14| into the PSG represen- 
tation of Eq. ( |A4| : GnMr, = mG T2 G Tl G c ^ G C ' 6 G T2 



r] 5 G T *Gc 6 , and Eq. ( |A5| : G a G Tl = Ti 6 G T2 G a , G a G T2 
tjtGtxGo-, we find 



r?! 1 X2 r]4Gc 6 (xi,X2) 



Gc 6 (xi + l,x 2 ) 
Gc 6 (xi,x 2 + 1) = rj^ Xl ri 5 Gc 6 (x 1 ,X2), 

G a (xi + 1,X 2 ) = l]^ 2 T] 6 G a (x 1 ,X2), 

G CT (xi,x 2 + 1) = Vi lr nGa(xi,x 2 ), 
whose solutions are 

/~<t \ xifxi — 1)/2— X1X2 xi aj 2/ ^< /n n\ 

Gce^i,^) = % ?? 4 % Gc 6 (0,0), 

G CT (xi,a;2) = r^ i:E2 77g 1 77y 2 G (T (0, 0). 



(A18) 



(A19) 



However, it worth mention that 774, 775, % and 777 are not 
independent!^. Because in their equations, either Gt x or 
Gt 2 only appears once, which means if we fix the inter- 
unit-cell gauge in a different way such that Gt x —> —Gt x 
or Gt 2 ~Gt 2 , the above four 77's will be affected. But 
this does not affect the mean field ansatz, as all the ansatz 
are given in the bilinear form which are invariant under 
this Z2 gauge transform. Therefore, make use of this I2 
gauge freedom, one can set two out of the four 77's to 
identity, say 774 = 775 = a . 



Substitute Eq. (A14), Eq. (A17|, and Eq. (A19) into 



the PSG representation of the rest of the definition re- 
lations, we find some constrains between the 7/'s. For 
example, from Eq. ( A6 ) one can obtain 



G <j{A)G a {B) = G a {B)G a {A) = (r) 6 r, 7 ) 1+ 2 Via- (A20) 

The left-hand-side is independent of (x,y), so must the 
right-hand-side be, therefore we must have 775777 = ctq, 



which means 7/5 = 7/7. Similarly from Eq. ( |A7[) and 
Eq. ( A8 1 we find 772 = 773 = 77s and from Eq. (|A9[) we 



find 775 = ?7 6 , so eventually 774 = 775 = rjs = rj 7 = <jq. 

Now all the G g (xi,x 2 ) has been reduced to G g (0, 0) 
with in a single unit cell, concluded as follows 



G Tl (xi,x 2 ) = ih 2 , 

Gt 2 (xi,X 2 ) = (7 , 

G r (xi,X2) = G T (0,0), 

G C6 (x 1 ,x2) = r 1 i l{xi ~ 1)/2 - XlX2 GcM0), 
G a {x 1 ,x 2 ) = r 1 * lX2 G a (0,0). 



(A21) 



The remaining task is to determine G7-(0, 0), Gc 6 (0,0) 
and G CT (0,0) from the following equations 

G T {A)JCG T {A)JC = G T {B)JCG T {B)JC = rj S , (A22) 

G r (B)ICGc 6 (A)Gr(A)IC = miG Ce (A), (A23) 

Gt{A)K,Gc 6 {B)G t {B)K, = miGceiB), (A24) 

G r (B)ICG a (A)Gr(A)IC = t^G^A), (A25) 

Gr{A)KG a {B)G T {B)lC = rnaG a (B), (A26) 

G a (A)G a (B) = G a (B)G a (A) = mo, (A27) 

(Gc^G^A)) 2 = (G Ce (A)G*(B)) 2 = mm , (A28) 
(G C6 (B)GcM)) 3 = (Gc 6 (A)G C6 (B)) 3 = 771779. (A29) 

The solution of the above equations leads to 144 algebraic 
PSG's which will be classified below. There are only two 
remaining SU(2) gauge freedom: the local gauge trans- 
form on site A or site B in the unit cell. 

We start from the solution of G7-. Let Gj- = a a + 
ia\0\ +ia 2 a 2 +ici3a3 be the most general form of a SU(2) 



matrix (with G 



Plug into the left-hand-side of 



Eq. (A22|, one finds 



-a 2 + <4)a 
2ia 2 (a 3 ai + a a 2 — ai<r 3 ). 



(A30) 



So if 778 = —(To j the solution is 02 = il } o-o = o-i = 
03 = 0, i.e. Gj- = ±ia 2 . Note that i<j 2 K. as a whole 
is SU(2) gauge invariant, thus the remaining gauge free- 
doms are preserved (even though er 2 seems to be a spe- 
cial direction). While if r] S = <tq, the solution is 02 = 0, 
a o + a i + a i = 1- O ne can choose Gr{A) = Gr{B) = cr . 
In this case, the SU(2) gauge freedoms on both sites are 
fixed. 

Class (I): 7 /8 = -cr . Then G r (A) = ia 2 , G r (B) = 
7771402, where 7714 = ±<7o is a new sign factor. Substitute 
into Eq. ( |A23|A24|A25|A26[ ), one finds 7711 = 7/12 = -7714, 
with no restriction on Gc 6 and G a . Fix the relative 
gauge between sites A and B by G a (A) = <tq, then from 
Eq. JA27J), G a (B) = rj w . Plug into Eq. rfAggl), 



(A31) 



According to the sign of 7717713, the class (I) is further 
divided into two subclasses. 
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Class (I)(A): 3717713 = ctq. Then the solution of 
Eg.([A3T|) reads G C M) = <?o, Gc 6 (B) = 7715 . Substi- 
tute into Eq. (A29), one finds 7715 = 771*79. 



The solutions in the class (I) (A) are summarized as 



(A32) 



G T {A) = ia 2 , 
G r (B) = ii]u<72, 
G Ce (A) = °o, 

Gc 6 (B) = 77i V9, 
G a {A) = CT , 
G a {B) — ?7l0 t 



which are controlled by 771, 779, 7710, 7714, providing 2 4 = 16 
PSG's. 

Class (I)(B): 7717713 = -a . Then from Eq. fA3l| ), 
the general solution of Gc 6 is a linear combination of 
iai, id2, 2C3. Because the global gauge freedom has not 
been fixed, so using this freedom, one can set Gc 6 {A) = 
ia\. Further assume Gc 6 {B) = 771779 (ifTi co s Q\ + 
(702 cos 02 + ica sin 2 ) sin0i), and plug into Eq. (A29l, 
one finds cos30i = l,sin30i — 0, whose solution is 
B\ = 0, ±27r/3, and there is no restriction on 02. 

The solutions in the class (I)(B) are summarized as 



G T (A) 


= «C 2 , 


G T (B) 


= «?7l4CT2, 


G a {A) 


= c , 


GM) 


= ??10, 


GcM) 


= i<?\, 


GcM) 


= ViVQ'i' 7 ! exp(i(720ie 



(A33) 



which are controlled by 771, 779, 7710, 7714, 0i, providing 2 4 x 
3 = 48 PSG's. Here 2 £ [0, 2tt) is a free angle. 

CTass fill: 7 7 8 = <t Then G r (A) = G r (B) = t7 . 
Therefore Eq. ( |A23|A24|A25|A26[ ) become 



KGcM)K = r hl GcMl 
!CGcM)£ = miGcM), 
JCGM)^ = m2G a (A) 1 

KG a (B)^ = VuG a (B)- 



The general solution of KG g K = G g is G g = e 



the general solution of JCG g JC = —G g 
According to the sign of 7711 and 7712 
further divided into four subclasses. 

Class (II)(Alhjhi = 7712 = ab- 
solution of Eq. ( A34 ) reads G Ce {A) = 
77 1 7 79 e^ e2 , G a (AJ~= e"* 8 *, G a (B) = 
acco rding to Eq. (|A27|), 4 = -0 3 . 

= e 21 
= (7 q . then 
<jq, thus 2 : 



(A34) 



ia29 while 

i(J 2 9 



Eq. (A28l, we obtain 



which implies e 2 ^2(0i+o 2 ) = 



is G g — ia 3 e 
, the class (II) is 

Then the general 

e^ e \ GcM) = 
= 77i e ICT294 . Then 
Substitute into 

as(fli-fla) = r?ir?13) 



reduced to e 



ia 2 {0 1 +9 2 ) 



Eq. (|A29j) can be 
= -0i. While 0! 



and 03 are related by 



h = 3 



7/17713 = er , 

t/2 ?7i77i3 = -co- 



(A35) 



The solutions in the class (II) (Al) are summarized as 



G T (A) 


= c , 


G T (B) 


= c , 


GcM) 




GcM) 


= ViV9e 


GM) 




GM) 


= T)lO e ~~ 



(A36) 



which are controlled by 771, 779, 7710, 7713, providing 2 4 = 16 
PSG's. Here 03 can be any angle, and 0i follows from 
Eq. ( [A35] ). 

Class (II)(B1): —7711 = 7712 = cr . The solution of 
G (j is the same as the class (II) (Al), however the general 
solution of G Cfi becomes: Gc (A ) = ia 3 e i,72ei , G C M) = 



—rjirigia^e 1 ' 72 2 . From Eq. (A28) one finds 7717713 = — uq. 
While from Eq. ( |A29| ), e zi,J2 ^ e ^ = 1, thus (0i - 2 ) = 
0,±2tt/3. 

The solutions in the class (II) (Bl) are summarized as 



G T (A) 


= c , 


G T {B) 


= c , 


GcM) 




GcM) 




GM) 




GM) 





(A37) 



which are controlled by 771,779,7710,(01 — 02), providing 
2 3 x 3 = 24 PSG's. Here 02 and 03 can be any angles, 
and (0i -0 2 ) = 0,±2tt/3. 

Class (II)(A2): — mi = —7712 = a . Then the gen- 



eral solution of Eq. ( A34 1 reads Gc 6 (A) 
GcM) = -r/iwcae 3 ^ GM) = «c 3 e' 



2c 3 e' 



1<J 2 V 3 



l(J 2 Ui 



JGa{B) 

Then according to Eq. (IA27), 3 



"77io7cr 3 e'" 

94. Substitute into Eq. (A28), then combining with 



Eq. (A29l, one finds 02 = — 0i, and 0i and 03 are related 
by Eq. (|A35|). 



The solutions in the class (II) (A2) are summarized as 



GM) 


= c , 


G T {B) 


= c , 


GcM) 




GcM) 


= -77i77 9 7;(73e lCT2ei 


GM) 




GM) 


= -f lw ^ a29 \ 



(A38) 



which are controlled by 771, 779, 7710, 7713, providing 2 4 = 16 
PSG's. Here 03 can be any angle, and 0i follows from 
Eq. ( |A35| . 

Class (II)(B2): 7711 = —7712 = cto- The solution of G a 
is the same as the class (II) (A2), however the general 
solution of G Ce beco mes G gJA) = e^ 201 , G C M) = 
rn?] 9 e ia292 . From Eq. (|A28[), it is found that 7717713 = 



12 



-cr . And Eq. \A29\ gives e 3la ^ 6l+0 ^ = 1, so (8 1 + 2 ) 
0,±2tt/3 



second term, and use 
The solutions in the class (II) (B2) are summarized as are two reasons 



Here we may simplify the expression by dropping the 

There 



j — ^2a=oXi Xj a a- 

First, consider the dihedral group 



G T (A) 


= Co, 


G r (B) 


= c , 


Gc 6 (A) 




Gc e (B) 




G„(A) 




G a (B) 





(A39) 



which are controlled by 771,779,7710,(6*1 + 9 2 ), providing 
2 3 x 3 = 24 PSG's. Here 9 2 and 6 3 can be any angles, 
and (61+62) = 0,±2tt/3. 

Now all the 144 algebraic PSG's has been classified. 
Given Eq. ( 13 1, one can check 77s = G7-/CG7-/C = — ctq and 
= (Gc 6 (A)G a (B)) 2 = (-cr Cf ,cr CT ) 2 = -er , which 
match the criterion of the class (I)(B). So the PSG of 
Kitaev spin liquid belongs to the class (I)(B) with 771 = 
= -V10 = -Vu = c and 6 1 = 2tt/3. 

Finally our classification is related to the previous 
wor y28| m ^ e following table. The number of PSG's in 
the class (II)(ii)(B)(/3) was miscounted in Ref. 28 as 24, 
which should be 8 instead. 

TABLE II: Relation between the new classification and the 
previous one. 



This work 


Ref.Hi 


(I) (A) 


(I) (A) 


00(B) 


(I)(B) 


(II)(A1) 


(II) (i) (A) + (II)(i)(B)(a) 


(II)(B1) 


(II)(h)(B)(a) 


(II) (A2) 


(II)(ii)(A) + (II)(ii)(B)(/3) 


(II)(B2) 


(H) 00 (B)(0) 



D 2 = {1, e l7vS i , e™ s i , e l7rS i}, which is a symmetry of the 
model Hamiltonian, and should not be broken in the spin 
liquid or Fermi liquid phase. Under these D 2 operations, 
Majorana fermions undergo sign changes, say for exam- 
nip pi^^i ■ — ^ v9 V"^ — ^ — V"^ V? — ^ — V? V? — ^ V 3 

which can be seen from the behavior of spin operators 
Si -Si, Sf -> -Sf, Sf -> Sf. Then any term that 
change the flavor of Majorana fermions acquires a mi- 
nus sign under at least one of the D 2 operations. So the 
D 2 symmetry preserves the flavor of Majorana fermions, 

and terms like XiXf an d *e a/37 xf x] are n °t allowed. Sec- 
ondly, in the time reversal broken phase like the super- 
conducting phase, the D 2 symmetry is broken. But in 
this case the bosons condense to a state described by 
Eq. ( 28 ) , which does not support any boson gauge cur- 
rent, i.e. Tr (a 3 B\a a Bj) = (a = 1, 2, 3). So the second 



term in Eq. ( B2 1 can not make a contribution to the mean 
field Hamiltonian in any case, and thus can be neglected 
for the sake of simplicity. Therefore the electron hopping 
term can be written as 



F * = ~lE E iXiXjTrasBj (-ia ) B 5 



(ij) a=0 

3 



(B3) 



(ij) a=0 v=\ 



For the Kitaev spin coupling term Hj, we first rewrite 
the spin operator to match Kitaev's convention by com- 
bining it with the neutral gauge charge K = 0, 



(B4) 



Appendix B: Mean Field Decomposition 

Rewrite the hopping term on a single bond in terms of 
Fi and B; matrices as 



E(4 



Cj a + h.c. 



TrasCjCj 



-Tr azBlFjFjB 



According to Eq. ( 10 1 

3 

2 



(Bl) 



^i = ^E^>+ 

Q=0 

a=l \ 0,7=1 J 



(B2) 



The single-occupancy projector (1 — 5) is appended to 
project out the holon gauge charge terms in Kf . Physi- 
cally (1 — 8) represents the probability that one electron 
actually appears on site so that the spin operator can 



make a effect. Substitute Eq. (B4| into Hj K 



Hj K = J K (l 



(ij) 



(B5) 



where a denotes the type of the bond (ij). 

Then by introducing the mean field parameters in 
Eq. ( 23 ) and following the standard slave boson mean 
field approach, Eq. (|B3|) and Eq. ( B5 1 can be decomposed 



to the mean field Hamiltonian Eq. ( 24 ) through Hubbard- 
Stratonovich transform, with additional Lagrangian mul- 
tipliers to enforce the doping and SU(2) gauge con- 
straints. 



13 



Appendix C: Case of Finite J 

Starting from the t-J-Jji model, H — H t + i?HK with 



H t = -tJ2 vJ ia c J<y V + h.c. - /i^c}^, 

(ij)a irr 

^hk = J ^ Si ■ Sj — Jk S°~Sj. 



(CI) 



Following the SU(2) slave boson theory, introducing the 

mean field parameters: ufj = (ixfxf)i w ij — {wufijj}, 
one obtains the mean field Hamiltonian 

Hmf = E ( E U &XiX! + E WrM^ + h.c.)) 

(ij) \a=0 v=l ) 

+e(e«-me^V 



\Z=1 



v=l 



(C2) 



where the fermion bond strength reads 



+ (1 - <5) 2 U° (J K <U - J H (1 - M) (C3) 



,8=1 / 

and the boson bond strength reads 



(C4) 



Q=0 



The index a denotes the type of bond (ij). The boson 
chemical potential /i is chosen such that — 

SN. The gauge potential a' is adjusted to ensure the 
gauge singlet condition (if?) = 0. 



The mean field phase diagram can be obtained by solv- 
ing the mean field Hamiltonian Hmf self-consistently. All 
phase diagrams contains SCi phase with Chern number 
+1 at small doping limit and SC2 phase at large doping, 



separated by the first order transition at S c . Tab. Ill 
list the values of 5 C for different settings of t and J. 
J/Jk = 1/8 corresponds to a = 0.8 according to the 
convention J = 1 — a, Jk = 2a. We conclude that small 
Heisenberg coupling will not affect the phase diagram 
much on the mean field level. 

On the type-3 bond, parameterize the mean field 
ansatz by u° = u , u\, = u?- = u bl = u a , wjj +w"j, — 



Then 



tw 



U° = - — + (1 - <5) 2 ((J K - JK - 2Ju b ). (C5) 
TABLE III: Kitaev spin liquid-FL transition point. 



t/J K 


Jh/ Jk 


S c 


10 





0.064 


1/8 


0.056 


5 





0.12 


1/8 


0.11 


2 





0.22 


1/8 


0.20 



The evolution of fermion mean field parameters with dop- 
ing at zero temperature is shown in Fig. [4] The first order 
transition between the Kitaev spin liquid and the Fermi 
liquid phases happens when = (at this point, the 
X° band becomes completely flat and can not gain more 
energy from the magnetic interaction). It is found that 
the mean field solution follows w = S and u b — £<5/(3Jk) 
at zero temperature, then the first order transition point 
6 C can be roughly estimated from the equation 



t6c 
2 



(1 - 5 c y (J K - J)u a 



2Jt5 c 



(C6) 



Considering the case of J = and large t, the transition 
point will be simply given by 5 C = 2u a J^/t 1 where the 
value of u a ~ 0.3 can be determined by the mean field 
solution. 
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